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1. Introduction and Overview 

Uncovering the elusive mechanism of electroweak symmetry breaking is one of 
the main goals of present and future accelerators. In these lectures I will concentrate 
on the most popular and simplest of all proposed mechanisms for that breaking. It 
makes use of a sector of fundamental scalar particles and the breaking of SU(2) x 
U(l) is achieved spontaneously by the vacuum expectation values (vevs) of some of 
these scalar fields, which have in general non-vanishing electroweak charges. 

A generic prediction of these models is the existence of physical scalar particles, 
the Higgs bosons, remnant of the electroweak breaking searched for in accelerator 
experiments. The aim of these lectures is to give an introduction to some theoretical 
guide available for that search in the form of Higgs mass bounds. The first section 
starts with some general statements that can be made with the only assumption that 
the Higgs sector is weakly interacting. Particularly relevant examples of models, 
like the pure Standard Model or its minimal Supersymmetric extensions are then 
put in a clearer perspective by contrast with the general case. 

The precise computation of an upper bound on the mass of the lightest Higgs 
boson in the Minimal Supersymmetric Standard Model is the topic of Section 3, 
while 4 and 5 are devoted to an equally precise computation of a lower bound on 
the Standard Model Higgs mass, from studies of the effective potential structure. 
These bounds lie in a mass region especially appealing for Higgs searches in the 
near future and are thus very relevant. Some implications that would follow from 
a Higgs discovery in such mass region are considered in Section 6. 
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2. Limits from spontaneously broken symmetries 



2.1 General Sum Rules 

Let {fa} be the set of Hermitian spinless fields in the theory with a potential 
V(fa) not necessarily polynomial, but invariant under some continuous symmetry 
G (global or local): 

fa^fa + e^fa, (2-1) 

where the 6> a, s are the generators of G, antisymmetric in our Hermitian basis. 

We are interested in the case of spontaneously broken G, so that the minimum 
of V(fa) occurs at fa = with (9 a v)i ^ for some of the en's. Furthermore, we 
will be interested in the scalar spectrum of states in the broken minimum V{. The 
first derivative of V with respect to the fields fa will be zero at the minimum by 
definition. The second derivative, 



Ml 



d 2 V 



dfadcpj 



gives the scalar mass matrix while higher order derivatives give scalar self-couplings 



d 3 V 



dfad^jdfa 



d 4 V 



ijkl 



d4>id4>jd4>kd4>i 



Our starting point is the invariance of the scalar potential under the transformation 
(2.1): 

V((j>i + e a 0%<l>j) = V((f>i) + e a ^9f k cf> k + 



d(pi 



(2.2) 



from which we obtain the identities 



0(pi 



(2.3) 



Taking derivatives of these functions W a with respect to the fa fields and evaluating 
them at the minimum, relations among different parameters of the potential can be 
obtained. The first derivatives give 



dW a 

dfa 



= 0=>0 = M? j (d a v) i , 



(2.4) 



which is just the familiar Goldstone's theorem: for every generator that does not 
annihilate the vacuum, M 2 has a zero eigenvalue. We note at this point that 
Vi can be written in general as a sum of a G-singlet Sj and a non singlet Vi. As 
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9 a v = 9 a v+9 a s = 9 a v we will drop the tilde in the following. The second derivatives 
of (2.3) give a relation between masses and cubic couplings: 

= 0^0 = e ljk {9-v) t + \M\ 9% k , (2.5) 

while the third derivatives relate quartic and cubic couplings: 

d 3 W a 

d<f>d<f>kd<f>l = = + + em&fi + euM, (2.6) 

(sum over repeated indices is assumed). Now, it is simple to eliminate cubic cou- 
plings between (2.5) and (2.6) to obtain a sum rule between masses and (adimen- 
sional) quartic couplings 

Mf j [(9 a 9 /3 v) i {9^9 s v) j + (9 a 9^v) i (9 /3 9 5 v) j + {9 a 9 5 v) t {9 l5 9~<v) J } 

= f ij ki(0 a v) i (0f i v) j (P'v) k (0 s v) l . 2 ' 7 

This is the central relation of this subsection. It implies that knowledge about the 
dimensionless scalar quartic couplings can be used to relate some scalar mass to the 
G breaking scale set by Vi, which is the only other dimensionful parameter entering 
(2.7). We will have occasion to see similar mass sum rules later on. Now it would 
be interesting to particularize (2.7) to the electroweak gauge group breaking and 
this example will provide a clearer picture of the relevance of (2.7). 

2.1 APPLICATION: G = SU(2) L x U(1) Y -> U(l) em 

There are three broken generators that we choose to be 9 a = {T 3 , T + , T_}. 
Corresponding to these generators there will appear three Goldstone bosons (to be 
eaten by the gauge bosons). Those are given by the (complex) vectors 

where N, N' are normalization constants related to the gauge boson masses by 
M z = gN/ cos 9 < w and M w = gN'/y/2. 

Conservation of electric charge implies there are only three non-trivial mass 
sum rules to be derived from (2.7) corresponding to the choices of {# Q , 6@, # 7 , 9 s } 
with charges adding up to zero. The scalar mass matrix will have block diagonal 
form breaking up in different submatrices for the differently charged scalar particles. 
Using a mass eigenstate basis, Mf- = MfSij the following mass sum rules are 
obtained _ 

3 

fzw = £{[Mf]%^ + m\ +) m\ 2 } , j (2 . 8) 

3 

f w = Y J {m? ) ?B 2 j + [M] ++) ] 2 | jDj f} 
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where f z = fijkiZiZjZkZi, fzw = fijklZiZjW k wf , f w = f ijkl WiW*w k Wi , are the quar- 
tic couplings of Goldstone bosons and 

*j = ^(Tfoj, Bj = ^[(T(T + 1) - Tl)v] r 

In (2.8), superindices in mass matrices indicate the charge of the states. 
Exercise. Obtain the corresponding expressions for Cj and Dj. Then prove the 
following relations 



E A i v i = E = L A i B i > <»• 
j j 



P 2 ' 



where G F is Fermi's constant, (v^Cf)" 1 = (246Gel/) 2 , and p = M W /(M Z cos8 w )- 
It is easy to transform these sum rules into useful mass inequalities. Consider 
the first equality in (2.8). If every (neutral) mass eigenvalue M- is substituted by 



the lowest one M^J n it follows that 

3[Mil] 2 5>?</z. 

j 

Furthermore, using the relations (2.9) and the Schwartz inequality, it follows 

(E,A,v,r i _ i 



E4> 



j 



2 ' 



Then one obtains 



[Mil] 2 < \fzv 2 < \f z {V2G F )-\ 



(note that possibly large singlet vacuum expectation values were dropped from the 
Ufc's). By similar methods the following inequalities can be obtained from (2.8): 

r M (0) I 2 < { — l^n* fw \ ( 2 10) 

We can see explicitly here that provided the quartic couplings / are weak one expects 
always a scalar with mass fixed by the electroweak scale. This holds irrespective of 
how complicated the Higgs sector is. In the case of the minimal Standard Model, 
with electroweak breaking described by a single Higgs doublet, the above inequalities 
are the same. In fact, the inequalities are saturated and one finds the familiar 
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relation between the mass and quartic Higgs coupling of the Higgs boson, M% = 
Xv 2 . In the general case there will always appear some scalar state whose mass 
cannot be made much bigger than the electroweak scale without making large some 
dimensionless coupling. It is natural to refer to that state as the 'true' Higgs boson. 
In general, the masses of other scalar states can instead be made heavy by choosing 
large values of the mass parameters in the lagrangian. 

2.2 Quartic Polynomial Potentials 

One remarkable property of the sum rules (2.7) obtained in the previous sub- 
section is that the only assumption on the form of the potential was the requirement 
of invariance under the action of G. Then, these sum rules can be obtained even 
if the potential is non-polynomial (e.g. in supergravity, or in potentials including 
radiative corrections). However, it will prove illustrative to consider a somewhat 
less general situation to gain some intuitive understanding on the relation between 
the spontaneous breaking of some continuous symmetry and the mass limits implied 
by it. In this subsection we will then concentrate in cases where the potential is 
indeed polynomial and at most quartic in the scalar fields. From the spontaneous 
breaking of the symmetry G of the potential we will obtain sum rules different to 
those derived in the previous subsection. 

Let us look more closely to the potential along the direction of the symmetry 
breaking. For this purpose use a mass eigenstate basis {</>a}- We can decompose 
every 4>a in a G-singlet part, sa, and a non-singlet part Ha- Fixing all singlet parts 
to their vevs the potential for the remaining /i^-fields is trivially invariant under 
G. In particular, this implies that no linear terms in the /^-fields will be allowed. 
Now consider the breaking direction in ZiA-space. Call {ua} a unit vector pointing 
in the direction of the breaking: ua = (^a) /v and h the normalized field along that 
direction: h = ^a n Ah>A- All the G-breaking is then given by the vev of the state 
h, while all other orthogonal states in /i^-space will have zero vev. The potential 
along h can then be written in the form 

V(h) = V(0) + \fh 2 + ^ah 3 + \\ h h\ (2.11) 

The value of h at the minimum, (h) = v, can be related to the parameters in the 
potential via the minimization conditions as usual 

^ = 0^ ^ + av + X h v 2 = 0. (2.12) 
oh 

The second derivative at the minimum can be written as 

9 2 V v -\ 9 2 V v -v o 2 2 r> o \ 2 

— — = > n A n B = V M A n A = n +2<rv + 3\ h v 

dh 2 d(f> A d(f>B 
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or, using (2.12) 

M A n \ = vv + 2X h v 2 . (2.13) 

A 

This sum rule does not seem to be particularly useful because the right hand side 
contains some unknown mass parameter a. However it is simple to see that the 
quantity av is always negative in the true minimum of the potential (2.11) [the 
degeneracy of the two minima of (2.11) at a = is lifted by the term ah 3 . The 
true minimum will then correspond to av 3 < 0]. As a result, av can be dropped in 
(2.13) to obtain the mass limit 

£ M\n\ < 2X h v 2 [M min ] 2 < 2X h v 2 . (2.14) 

A 

The last mass inequality follows the same line of Sect. 2.1 and uses Y2 n A = 1- 
Note that now, the quartic coupling A^, is not necessarily related to the Goldstone 
couplings that appear in the mass inequalities of (2.10) implying that in general, 
the mass bound (2.14) is different from those. 

We can improve the previous derivation if the symmetry group G contains an 
SU (2) subgroup (this apply in particular to the electroweak gauge group and is the 
reason why we concentrate on it. One can construct a similar line of derivation for 
any U(l) subgroup of G). In that case, the Hermitian scalar fields 4>j will belong 
to some SU(2) multiplet of dimension 27} + 1. We can define the following parity 
transformation 

P: •( l) 27 ^. 

This transformation can indeed be defined for all the fields in the theory and also for 
products of any number of fields. Call even the fields (or product of fields) invariant 
under P and odd the fields that change sign under P. In particular, SU(2) singlets, 
triplets, etc are even while doublets, quadruplets, etc are odd. 

By fixing all even fields to their vevs we can obtain the potential for odd fields 
only. In the odd-field space we then define the normalized field 

° odd 

along the direction of the breaking. Again, all combinations of odd fields orthogonal 
to <f have zero vev and the potential V(<p) must be P-invariant. Then, not only 
linear terms in <p are forbidden, but also cubic ones: 

V(<p) = V(0) - imV + 

If tp were a mass eigenvalue it would have mass M 2 = X^v 2 . In the general case 
it is clear that some eigenvalue will have mass below that value and the following 
inequality follows: 

[M min ] 2 < X^v 2 . (2.15) 



7 



Note that, in general, both the quartic coupling and the vev appearing in this mass 
bound will be different from the ones derived previously. 



2.2 APPLICATION: Lightest Higgs in SUSY extended models 



It is straightforward to apply the mass bound (2.15) to derive a tree level upper 
bound on the mass of the lightest neutral Higgs boson in general supersymmetric 
models. The quartic coupling for <p is obtained from two sources: F and D terms. 
The only terms in the superpotential than can contribute to A v through F-terms 
will have the form W3 ~ /i$ e $ $ (the subindices indicate whether the chiral 
superfields $ are even or odd under P). The contribution of such terms to V(ip) 
will be then 



8 F V(<p) 



2 2 4 



(2.16) 



where the dimensionless Yukawa couplings h^s are in general accompanied by some 
rotation angles coming from the projection of the odd fields into the ip direction: 
<k = Vi/v a . The contribution of D-terms depends on the gauge group. Assuming 
the minimal SU(2) L x U(l)y it is easy to obtain 



SnV(<p) = -(g 2 + g' 2 ) 



i,odd 



(2.17) 



Adding (2.16) and (2.17), eq. (2.15) gives the following mass limit 

2 



M^<(g 2 +g") 



i,odd 



P + 2j2hMvl 



(2.18) 



In softly broken SUSY models, the details of the supersymmetry breaking do 
not affect the quartic couplings (at tree level) so that the scale of the mass bound 
(2.15) will be fixed only by the Fermi constant. As we have seen, this is not a 
particular feature of supersymmetric models but holds with complete generality. 
The Minimal Supersymmetric Standard Model has the particularity that only gauge 
couplings appear in (2.18). As these are measured experimentally, a well defined 
tree level mass bound results. In general, numerical values for the bound (2.18) 
can be obtained only by setting upper limits on the (unknown) hi couplings, e.g. 
by triviality arguments, requiring that they remain perturbative up to some large 
energy scale. 

Exercise. For spontaneous breaking of the electroweak gauge group: a) if driven 
only by 577(2) doublets, show that the bounds (2.10,2.14,2.15) coincide; b) if the 
scalar potential has a SU (2) custodial symmetry after the breaking, show that the 
bounds (2.10) are equal. 
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2.3 Interplay between different bounds and the decoupling limit 

All mass bounds derived so far have in common the following property: they 
follow from a sum rule of the form (a simply counts different bounds) 

(ifaW 2 ^) = X a vl (2.19) 

where <p a and 4> a are certain scalar fields [that can in general be normalized such 
that ((p a \(f>a) = 1] ) M 2 is the scalar mass matrix, A a is some scalar quartic coupling 
and v a is of the order of the breaking scale. We will call (2.19) the 'vector-form' 
of the bounds. From them, the usual 'scalar-form' of the mass bounds are de- 
rived straightforwardly by repeated use of I = YIa l^-X^-l' where the \A) are mass 
eigenvalues: 

(¥>«|M 2 |0 O ) = X a v 2 a = J2MA)M 2 (A\cf> a ) > M 2 Y,MA)(A\$ a ) = M 2 . 

A A 

In table 1 we present all five mass bounds obtained previously in vector and scalar 
forms. The first column gives the vectors entering (2.19). When only one vector is 
shown it is assumed <p = 0. Note that only the second bound has different vectors. 
Exercise. If ((p\M 2 \(p) < X^v 2 and {(j)\M 2 \4>) < A^J with \<p) ^ k\4>), then 
{(f\M 2 \(j)) has also a bound of the form Xv 2 . a) find A and v 2 ; b) apply this to 
(2.10). 



I<M, \<Pa) 


K 




M 2 h < 


<f>z ~ E 4 tliMi 


2^Aoo(G ) 4 


E^ 2 


\Xq V 2 




^Xq c \g°g + 2 


£—G~ x 


Xq c v 2 




\X CC \G + 4 




\X cc v 2 


I Ei Mi 


|A0 4 


£,«< 2 


Xv 2 


~ Q Ei,odd V i ( t ) i 


|A ^ 4 


V v 2 


Kv 2 



Table 1: Summary of mass bounds. 



In the following we will always assume the equality <p a = (j) a . The second 
column gives the relevant quartic couplings by showing the corresponding piece of 
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the scalar potential. G°, G ± are the Goldstone bosons. The third column gives the 
values of the vev's Vi and the fourth lists the scalar form of the mass bounds. 

The vector form of the bounds contains extra useful information not contained 
in the scalar form. As we have shown, in a particular model one can obtain several 
different mass bounds for the lightest Higgs scalar and by choosing the stronger 
among those one gets the best limit. Here we will show how the vector form can be 
used to extract further information on the scalar spectrum by the interplay between 
different bounds. 

Consider the inequality 

J2m)Mi(A\0) < Xv 2 . 

A 

If some state \A) is much heavier than the breaking scale, M\ 3> Xv 2 , its overlap 
with \4>) has to be correspondingly small: 

MA)\ 2 <^ 2 /M 2 A ^Q. 

This can be the case if all mass parameters in the potential are made very 
heavy (~ M) while the pattern of symmetry breaking is held fixed (vi constant. 
The required fine-tuning is bigger the heavier the mass scale M is made). We will 
call such situation the decoupling limit. In such limit one expects that all scalars 
will have masses of the order of the heavy scale M 3> v with the exception of the 
state constrained by the mass limit Xv 2 (and possibly others whose mass is protected 
by some symmetry). Then, if all scalars but one are very heavy, the light state |1) 
will satisfy 

Ml)\ 2 = l-0(v 2 /M 2 A )^l. 

This means that the state \4>) appearing in the vector- form of the bound is precisely 
the state that will remain light in the decoupling limit. 

Exercise. Given that |1) — > \<f>) in the decoupling limit, can you prove that M 2 — > 

Xv 2 = (0|M 2 |0)? 

Suppose now that your model has two 'linearly-independent' bounds Aii> 2 = 
(^i|M 2 |(/>i) and X 2 v\ = (0 2 |M 2 |0 2 ) with ^ n\(j) 2 ). What would be the lightest 
state in the decoupling limit, \4>\) or |<^> 2 )? Of course the way out of this paradox 
is that the pure decoupling limit cannot be realized: at least two states will remain 
light: 

Exercise. In the situation described above, prove that the mass of the second-to- 
lightest state M 2 satisfies the inequality 

M 2 - M 2 < — ^- { JXxvl-Ml + Jx 2 v 2 -Ml 
sin a y v v 

where a is the angle between |0i) and \<fi 2 ). 
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The existence of a bound on the mass of the second to lightest scalar in such a 
situation can be easily understood in terms of the following geometrical construc- 
tion. Consider the mass ellipsoid in the (multidimensional) scalar field space. Its 
axis lie along the directions of the eigenvalues with length proportional to the inverse 
of the corresponding squared mass eigenvalues. One mass bound of the form we have 
considered is given in this picture by some vector with length (At> 2 ) -1 and direction 
determined by The fact that the bound is satisfied implies that the ellipsoid 
can at most touch the tip of this vector but cannot intersect it: (<p\M 2 \(p) < Xv 2 . 
This is represented in the 2D example of fig. 1 a. As is clear form the picture there 
is necessarily an eigenvalue with mass M| < Xv 2 . By choosing all mass parameters 
heavy while keeping v fixed we can flatten the ellipsoid (i.e. make the mass eigenval- 
ues heavy) except in the vector direction of the mass bound, where the flattening is 
obstructed by the presence of the vector bound, that the ellipsoid cannot intersect. 
This is shown in fig. 1 b. From it, is clear that when all mass eigenvalues but one 
are heavy, the light state is given by \<f>). Finally, in the presence of two linearly 
independent mass bounds, the ellipsoid can be flattened at most to an ellipse (see 
fig. 1 c). The two light states, undetermined in general, will be linear combinations 
of \<pi) and |02)- Obviously, this can be generalized to any number of independent 
mass bounds and one will expect always a number of light states in correspondence 
with the dimensionality of the vector space {|0 a )}- 




Fig. 1: Geometrical representation of the working of mass bounds. 

Exercise. Suppose that there is some 'triplet impurity' in SU(2)l x U(1)y break- 
ing, i.e. some 577(2) triplet takes a vev x and contributes to the weak vector boson 
masses in addition to the usual doublets. The vev x is bounded to be small by the 
Ap constraint. In the decoupling limit, what triplet admixture do you expect in 
the light Higgs boson? Of which order would be the second to lightest Higgs mass? 

2.2 Production Cross Section of the light Higgs boson 

Concerning electroweak symmetry breaking, besides the knowledge about lim- 
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its on light Higgs masses one would be interested in determining whether such 
states can be detected experimentally at all. It is here that the information on the 
composition of the light Higgs becomes crucial because it determines directly the 
production and decay cross sections. Here we will limit ourselves as an example, 
to the Higgs-strahlung production mechanism of a light Higgs boson in a e + e~ 
collider: e + e~ — > Z* — > Zh±. The cross section for this process is proportional 
to the corresponding cross section for a Standard Model Higgs boson of the same 
mass. The proportionality coefficient is fixed by the gauge properties of hi. More 
precisely, the coefficient is given by the overlap between hi and the \<j>z) state listed 
in table 1: 

, . . l2 a(e + e~ — > Zhi) 
a(e+e -> Zh S M) 

Suppose that the model contains a scalar Higgs which is a gauge singlet. If the 
light state hi has a large overlapping with the singlet the production cross section 
for hi will be reduced. In this situation even if hi is forced to be below some mass 
bound it may be very difficult to produce and detect in accelerators. However, the 
important point is that from the vector form of the first bound listed in table 1 some 
state with a no n- vanishing overlap with <pz must remain light. In other words, if 
the lightest Higgs turns out to be orthogonal to 4>z one can still use the information 
from (<pz\M 2 \(pz) to put a bound on some other Higgs. Technically this is realized 
in the following way: consider the quantities (4>z\M 2 — m 2 N \(f) z) , where N numbers 
the (scalar) Higgs mass eigenvalues, mi < mi-... For N = 1 one obtains 

^z\M 2 \$ z ) = Xv 2 < 0, 
which is the original bound. For N = 2 the following inequality results 

m 2 2 < 

From this relation we see that if hi becomes singlet dominated and 

m\ < Xv 2 , i.e. the second to lightest eigenvalue satisfies the original bound. If on 
the other hand hi — > <fiz then no bound on ?ri2 can be set. In general, for the N th 
eigenvalue one finds 

2 Xv 2 - m 2 S 2 N 

m N - \ ' 
N 

with 

JV-l 

p =i 

When Sn is small, for the first N — 1 light Higgses having reduced couplings to the 
Z, the bound on m 2 N is then stronger. This effect can ensure sufficient production 
of some light scalar provided there are not too many singlets. 



Xv 2 - m 2 \{(j)z\, 


hi)? 


1 - 


(<t>z\ 


hi) 


2 
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3. The lightest Higgs boson in the MSSM 

3.1 The MSSM Higgs sector at tree level 

Supersymmetry requires that the minimal supersymmetric extension of the 
Standard Model, MSSM, contains two Higgs doublets (of opposite hypercharge) to 
give masses to all quarks and charged leptons. The most general tree level potential 
for the Hi, H 2 Higgs doublets, gauge invariant and renormalizable is then: 



V = m\\Hi\ 2 + m 2 2 \H 2 \ 2 + [m\ 2 Hi ■ H 2 + h.c] 



+ -^^2\H 2 \ 



+ A3 1 Hi 1 2 1 H 2 1 2 + A4 1 Hi ■ H 2 



(3.1) 



+ 



1 



A 5 (#i • H 2 f + X 6 \Hi\ 2 (Hi ■ H 2 ) + X 7 \H 2 \ 2 (Hi ■ H 2 ) + h.c. 



with 




Hi 
H° 2 



The quartic couplings in this potential are constrained by supersymmetry to be 

-if- A 5 = A 6 = A 7 = 0, (3.2) 



\i = \2 = \(9 2 +g% Xs = \(9 2 -9% A 4 = -lr 2 



i.e. they are given in terms of the gauge coupling constants. The projection of (3.1) 
on the neutral Higgs components gives the potential 



V(H°, tf 2 °) = mj\H°\ 2 + m\\Hl\ 2 + [m\ 2 H^H% + h.c] 
+ 1 -(g 2 + g' 2 )(\H 1 \ 2 -\H° 2 \ 2 ) 2 . 



(3.3) 



The minimum of this potential determines the vevs (H®) = vi/y/2 and (H 2 ) = 
V2/V2 with v\ + v 2 = (246 GeV) 2 fixed by the gauge boson masses. The ratio 
tan (3 = v 2 /vi is a free parameter. The Higgs spectrum in the broken minimum just 
described consists of two (CP even) scalars h°,H°, one (CP odd) pseudoscalar A 
and two charged Higgses if ± . Two of the three mass parameters in (3.3) can be 
traded by vi and v 2 , so that, at tree level, the properties of the Higgs sector (masses, 
mixing angles and couplings) are determined by one mass parameter (usually taken 
to be the mass of the pseudoscalar, uia) and tan/?. 

The discussion of Section 2 should have made clear that relations (3.2) have 
direct and important consequences for the Higgs spectrum. The tree level masses 
for Higgs bosons are: 



m 



H ± = m w +m A , 



m 



h,H 



m 2 A + m| =F \j (m 2 A + m|) 2 - Am 2 z m 2 A cos 2 2/3 



(3.4) 
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The two neutral Higgses are linear combinations of H® r 2 = y/2ReH® 2 , see fig. 2 




Fig. 2: Electroweak symmetry breaking in the MSSM. 

H° = (H° r - Vl ) cos a + (if£ r - u 2 ) sin a, 
/i° = (H% r - v 2 ) cos a - (H° r - Vl ) sin a, 

where the mixing angle a is given by 



cos 2a = — cos 2(3 



™\-m* z 

2 2 ' 

m H~ m h 



sin 2a = — sin 2/3 



m 2 H + m\ 
m 2 H -m 2 h 



The couplings of H° and h° relative to the Standard Model Higgs boson are then 
determined by a and (3 as shown in table 2. 





WW,ZZ 


ZA 


tt 


66, r + r~ 


h° 


sin(/3 — a) 


cos(/3 — a) 


cos a/ sin /3 


— sin a/ cos /5 




cos(/3 — a) 


sin(/3 — a) 


sin a/ sin/? 


cos a/ cos /3 



Table 2: Neutral MSSM Higgs couplings relative to the corresponding Standard 
Model Higgs coupling: gsuSY /gSM- 



We can apply some of the general results examined in section 2 to the MSSM 
case just described. For instance, if the only mass parameter available in the Higgs 
sector, rriA, is made heavy (compared with the electroweak scale as given by Mz), 
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eqs. (3.4) show that all Higgses will acquire masses of order tua except one neutral 
Higgs that remains light, with 



m 2 h < Mf cos 2 2(3. 

The inequality holds for any value of rriA and is saturated when rriA — > oo. This 
bound explicitly derived here is the particularization to the MSSM of the general 
bound (2.18) and is associated with a vector \<p) in (H°,h°) space that lies along 
the direction of the breaking, see fig. 2. In the decoupling limit, tua n%z, we 
know that the light state h° will tend to be aligned with the vector \ip). That is, 
a ^ (3 — 7r/2or 

h° ~ (ui/i? + v 2 h° 2 )/v = h\ cos /3 + h% sin/?. (3.5) 

In addition, a look at table 2 shows that the couplings of the light state h tend to 
the Standard Model values in this decoupling limit. Such a Higgs scalar should be 
detectable then at LEPII (when uia ~ mz the coupling h?ZZ can be suppressed 
closing the Higgs-strahlung channel but then the complimentary production mech- 
anism e + e~ — > Z* — > h°A° becomes important). 

In summary, a tree level analysis of the MSSM Higgs sector predicts the ex- 
istence of a CP even Higgs scalar with mass below Mz that should be detectable 
at LEPII. This would represent a stringent test for the simplest realization of the 
supersymmetric extension of the Standard Model. However, as is well known, the 
effect of radiative corrections modifies this expectation in a dramatic way. 

3.2 Radiatively corrected m^o . Dominant effect 

After including radiative corrections, the mass of the MSSM lightest Higgs 
boson is no longer determined only by rriA and tan (3 but will depend on the rest of 
parameters of the theory. In particular, the most important corrections come from 
top-stop loops that shift the squared mass of h° by 

Here M t and are the top and stop masses respectively. Some comments are in 
order. Note first that the contribution cancels if top and stops were degenerate as 
would correspond to the supersymmetric limit. Also note the strong dependence 
with Mt. As an example, if Mt = 170 GeV, Mf = 1 TeV, the maximum tree level 
value m h o = Mz would be shifted by (3.6) by as much as 30 GeV (the increase can 
be even more dramatic when the tree level mass is smaller). 

The radiatively corrected Higgs sector of the MSSM has been intensively stud- 
ied during the last years using three main tools: diagrammatic calculations, effective 
potential techniques and renormalization group methods. It is instructive to obtain 
(3.6) using these three approaches. 
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Exercise. Calculate the dominant correction (3.6) diagrammatically. Assume for 
simplicity rriA 3> w^o, i.e. work only with h as defined in (3.5). Consider the 
diagrams of fig. 3 using the rules h°tt: H t /V2; h°h°t* L R t L ,R- H h where H t = 
h t sin/3. Also, M 2 = H?v 2 /2, M? = m 2 + M 2 + 0{g 2 v 2 ), M? = ml + M 2 + 
0(g 2 v 2 ). [Left-right mixing in the stop sector can be neglected to compute the 
dominant effect (3.6)]. 



h -o- ----- -o ' 

Fig. 3: Top and stop one-loop diagrams contributing to m h0 . 
Exercise. The one-loop MSSM effective potential when tua 3> mz is given by 

V(h) = -\m 2 h 2 + \\h A 
2 8 



6 

+ 



64tt 2 



M? q q M 2 3 ' 

<dog^ - §) + M*(.og-|, - §) - 2M*( log ^ - |) 



with A = [(g 2 + g' 2 ) /A] cos 2 2/3 and M? ^ = if left-right mixing is neglected. 

The second derivative of the potential in the minimum ((f)) = v will give the one- loop 
corrected Higgs mass. 

The third method, perhaps the most elegant, is based on an effective field 
theory approach: when the SUSY spectrum is around some scale Ms 3> Mz (in 
particular Mf ~ tua ~ Ms) the effective theory below Ms is just the Standard 
Model. In particular remember that ttia 3> M z implies that the light Higgs boson 
[given by (3.5)] has Standard Model-like properties. The quartic coupling of h° is 
determined by supersymmetry at the scale Ms: 

X(M S ) = ^[g 2 (M s )+g' 2 (M s )]cos 2 2(3, (3.7) 

and its value at the electroweak scale can be computed integrating the renormaliza- 
tion group equations in the effective theory between Ms and the low energy scale 
(e.g. M t ) with (3.7) as a boundary condition. The mass of h° can then be obtained 
from A at the weak scale by using the Standard Model relation m\ = Xv 2 . It is 
then straightforward to reinterpret the result (3.6): the logarithmic dependence 
arises from the running of A from Mj ~ Ms to Mt and the M t 4 dependence comes 
from the dominant piece of the Standard Model beta function for A: 

P SM 12 4 
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Note that this renormalization group method is particularly well suited to study 
the upper bound on m? which is obtained in the large uia limit. 

A few comments on the structure of the radiative corrections to m? are in 
order. The fact that one-loop radiative corrections to this mass can be very sizeable 
(and even larger than the tree level mass if tan/3 is small) does not mean that 
perturbation theory is not reliable. A large ratio of one-loop corrections to tree- 
level contributions arises because the tree-level result does not depend on the large 
top Yukawa coupling while one-loop corrections do. Furthermore the one-loop result 
is enhanced by a logarithm of a possibly heavy mass (Ms) to a light mass (M t ). The 
loop contributions in the adimensional ratio Am^o/M 2 are basically of the form 



k a n ~\ 



n k=0 



where a = hi /An 2 and log ~ log(Mj/M t 2 ). The terms for a given n come from 
the n-loop corrections, with the index k corresponding to hierarchically organized 
contributions (provided the log is sizeable): leading-log terms (k = n), next-to- 
leading (k = n — 1) and so on. For the perturbative expansion to be reliable we 
should require a log < 1 which is satisfied for the current top mass values. The 
use of renormalization group techniques permits to reorganize the loop expansion 
resumming to all loops the numerically most important corrections (leading, next- 
to-leading, etc). In the next subsection we discuss how to implement a computation 
of the loop corrected mass m^o that will include up to next-to-leading radiative 
corrections. Use will be made of the three methods just sketched above for the 
dominant correction. 

3.3 Radiatively corrected m^o . Next-to-leading log computation 

The ingredients for a next-to-leading log computation of the radiatively cor- 
rected will be discussed in this subsection. We assume that the supersymmetric 
spectrum can be described by a common mass Ms well above the electroweak scale. 
In particular vtia ~ Mf ~ Ms and below Ms the effective theory is the Standard 
Model. As we saw, the quartic Higgs coupling at the weak scale will determine the 
light Higgs mass, while at the supersymmetric scale its value is fixed by supersym- 
metric parameters. First, we can see that the integration of the coupling A from 
Ms down to the electroweak scale indeed resums some series of corrections to all 
loops. For example, from dX/dt — (3\, where t = log/U, we get 

j- M s 

X(M s )-X(M t )= / (3 x (t)dt 



M S 



ji=M t 



d(3 x 



P\(t ) + (t- to)^r(to) + - + -r(* - *o) 



d n (3 } 



(3x (t ) log 



M s , I dp, 



+ 



2 dt 



dt 



(to) 



log 



Ms_ 

M t 



+ ... + 



log 



dt n 
M s 



+ 



dt. 



Mt 



+ 
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Inserting in the above expression a loop expansion for the beta function one recovers 
the general estructure of the radiative corrections discussed at the end of the previ- 
ous subsection. The one-loop approximation for the beta functions corresponds to 
the (all loop) leading log contributions. Using two loop beta functions would resum 
also the next-to-leading logs and so on. 

The fact that information on higher loop corrections can be obtained with the 
knowledge of just one-loop beta functions is due to the magic of the Renormalization 
Group. Let us have a closer look to it from the point of view of the effective potential 
V(ip). The starting point is the observation that a change in the renormalization 
scale — > ji + dji doesn't change the physics. The invariance of the potential under 
such change V({p((i), A^/u); y) = V((p(/i + d/i), + dfi); [i + dfx), can be expressed 
in differential form as (sum over repeated indices implied) 



' d d d 

o\i dip d[i 



V = 0. (3.8) 



Here A^ stands for a generic coupling with corresponding beta function fli, <p is the 
Higgs field, with anomalous dimension 7 (792 = —d(p/dt). Note that beyond tree 
level the effective potential also depends explicitly on \x through logarithms. Then 
(3.8) connects contributions to the potential of different orders in h. For example, 
at order h 1 one has 

„(i) dV (1) dV dVi 

where Vq and V\ are the tree-level and one loop potentials respectively and the 
index (1) in the (3 and 7 functions indicate one-loop approximations. Eq. (3.9) 
tells that knowledge of Vq, and one- loop rg functions allows the computation of the 
leading-log one- loop contribution in V± (which goes like log n) . 
Exercise. The reciprocal is also true. From the Standard Model effective potential 



% 



/I 2 



with n t = -12, Ml = h 2 ip 2 /2; n w = 6, = g 2 (p 2 /4; n z = 3, M| = (g 2 + 

g' 2 )¥ 2 /4; n h = 1, M 2 = 3\cp 2 /2 - m 2 ; n x = 3, M 2 = Xcp 2 /2 - m 2 (Goldstone 

bosons), obtain 0$ and ^ knowing that 7^) = S[h 2 - (l/4)g' 2 - (3/4)^ 2 ]. 



Defining the operators 



p(„) =Q(n)d__ („) d 



we can write the fi 2 expression of (3.8) as 



v^v + v^v 1 + ^ = o, 
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which would imply that knowledge of rg functions to two-loop order permits to 
obtain the leading and next-to-leading two-loop contributions in V<i (provided Vo 
and Vi are also known). 

The procedure can be extended to order h n : 

V (n) VQ + V (n-D Vl + ... + p(l)y n _! + V— ™ = 0, 

and recursion relations can be written for the nth-loop leading V^ L and next-to- 
leading log V^ TLL contributions: 

V (2) V LL + V (l)yNTLL ?K_ = Q 

n 2 n-l 

From them we see that V£ L can be obtained from Vq and T>^> while to obtain 
V^ TLL one needs in addition Vi and V^ 2 \ The general statement is that the L th 
loop potential, with parameters running with L + 1 rg functions resums contribu- 
tions up to L^-to-leading order. In our particular case we shall use the one-loop 
effective potential with parameters running at two-loops to resum leading and next- 
to-leading corrections. This should be compared with the approximmation in which 
the tree level potential with parameters running at one- loop is used, which would 
resum only the leading logs. 

An approximation for the effective potential truncated at some loop order, like 
the one we are using, will have some residual scale dependence. Understandably, 
this dependence will be ameliorated if we use a one-loop expression for V with two- 
loop running parameters as compared with the tree-level (one-loop rg improved) 
potential approximation. The goal is to compute reliably the second derivative of 
the potential (to be related with the Higgs mass) in the electroweak minimum, so 
that one would like to keep control on the scale dependence in that region of the 
minimum. 

One way of doing this is the following. If the potential were known exactly it 
would be exactly scale independent. In such a case it is simple to show that the vev 
(<fi(t)) should run with the scale in the same way as the field <p(t) does 



<p(t) = y? c exp io 7(t )dt = (p c £(t), 



where <p c is the classical field. Then, the ratio (<p(t))/£(t) gives a measure of the 
scale independence of the potential used to extract (ip(t)): the ratio should be flat 
in the region where the potential is more scale independent. This fact can be used 
to determine numerically some scale t* where indeed the above ratio turns out to 
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be flat. This is shown in fig. 4 , where the ratio (<p(t))/£(t) is plotted as a function 
of the scale for two different approximations for the effective potential: the dashed 
line corresponds to the tree-level (one-loop rg improved) effective potential and the 
solid line to the one-loop (two-loop rg improved) case. The figure shows clearly the 
improvement in scale independence if the second approximation is used. The scale 
where the corresponding curve is flatter determines the scale t* : not surprisingly it 
is of the order of the mass scale involved in the problem i.e. the top mass. 

As a reassuring cross check one can also compare the real running of second 
derivatives (at the minimum) of the potential used with the running if the potential 
were scale invariant. Again it is simple to prove that d n V/d(p(t) n runs like £~ n (t) 



in the latter case. Then the ratio of ml 



h,der 



d 2 V/dip(t) 2 at (<p(t)) to ml(t) 



m h(t*)£ 2 (t*) / £ 2 (t) measures the scale independence of the potential approximation. 
It turns out that this ratio becomes flatter at approximately the same scale t* 
previously found which is then used as the best and most reliable scale choice for 
numerical computations. 
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Fig. 4: Ratio (<£>(£))/£(£) as a measure of scale independence for the effective 
potential in two different approximations as explained in the text. 

Once the scale dependence of the potential is taken care of, one can compute 
the second derivative at the minimum and relate it to the Higgs mass. This mass 
will then include logarithmic corrections up to next-to-leading order. However one 
should realize that this is not yet the physical mass, i.e. it does not correspond 
to the pole of the Higgs propagator, Tn(p 2 ) = p 2 — (m 2 R + Ur(p 2 ), where p is the 
external momentum, mjt is the renormalized mass and the one-loop self-energy. 
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Remembering that the effective potential generates 1PI diagrams with zero external 
momentum we get 



m 2 h (t) = 



d 2 V 



-T R (p 2 = 0)=m 2 R + U R (0). 



The pole mass being defined by T R (p 2 = Mjj) = 0, we arrive at 
M 2 H = m 2 (t) + U R (p 2 = M 2 H ) - U R (p 2 = 0). 



(3.10) 



It can be shown that the scale dependence of m|(t) cancels at one- loop with the 
self-energy difference correction giving rise to a pole mass scale independent up 
to higher orders. The diagrammatic calculation of the one-loop self-energies adds 
to the Higgs mass the one-loop corrections not accesible to rg resummation. A 
similar self-energy correction should be included to relate the running top mass 
mt{t) — ht(t)v £(£)/\/2 with the top pole mass M t . The dominant piece comes from 
QCD radiative corrections. In MS it reads 



M t = m t (fi = M t ) 



1 4 a,(M t ) 

3 7T 



(3.11) 



The last piece for a consistent computation of m^o at next-to-leading order 
corresponds to the inclusion of one-loop threshold corrections to the boundary con- 
dition (3.7). They arise when integrating out the heavy supersymmetric particles. 
The dominant contribution corresponds to stops and is represented diagrammat- 
ically in fig. 5 . It is a simple Exercise to compute it either by expanding the 
contribution of stops to the MSSM potential in powers of the background Higgs 
field or by diagrammatic calculation. The correction to A(M$) is proportional to 
the stop mixing M? ~ = H t <pX t , where H t = h t sin (3 and X t = (A t + fx cot (5): 



dX -^Ml 



1 - 



X 



2 1 



12MJ 



(3.12) 
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Fig. 5: Two types of supersymmetric diagrams giving the dominant threshold 
corrections to the quartic Higgs coupling. 
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The shift in A reaches a maximum value for X? = 6M| that corresponds then 
to the maximum of the Higgs mass (maximal mixing case). The case of negligible 
mixing X t ~ (minimal mixing case) will in general correspond to the minimum 
value of the Higgs mass (when the rest of parameters is fixed). Note that the 
correction (3.12) can be negative if X? > 12MJ. However in that region of param- 
eters is easy to run into problems with color or charge breaking minima in the full 
supersymmetric scalar potential. 

This completes the list of ingredients for the next-to-leading log computation 
of m/jo. An example of the results is plotted in fig. 6 which gives the (physical) 
Higgs mass versus the (pole) top mass. 
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Fig. 6: Upper limits for the lightest Higgs boson in the MSSM with the scale of 
supersymmetry Ms = 1 TeV . Solid lines: tan/3 ^> 1, dashed lines: tan (3 = 1. 
For a given tan (3 the upper (lower) curve corresponds to Xf = 6M J (Xt = 0) . 

The scale of the supersymmetric particles is set to 1 TeV which is roughly the 
upper limit from naturality arguments. As the pseudoscalar mass via is then much 
larger than the electroweak scale, the masses shown in the figure are actually the 
upper limits for m^o (for the corresponding value of Ms- The masses do increase 
logarithmically with Ms)- Solid lines correspond to the case of large tan (3 while 
dashed lines have tan (3 = 1. In each pair, the upper curve is the one for maximal 
mixing X^ = 6Mj (giving then the absolute bound on m^o) and the lower for the 
minimal mixing case X t ~ 0. Short-dashed vertical lines give the CDF/DO range for 
the top mass. From the figure is clear that LEPII can miss the lightest Higgs boson 
if e.g. tan (3 turns out to be large and stops are heavy with substantial mixing. 
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A few comments are in order before closing this section. When there is still a 
hierarchy between Ms and M t but uia ~ M t one can repeat the above procedure 
with the difference that now, the effective theory below Ms is a two Higgs doublet 
extension of the Standard Model. The Higgs potential is then that written in (3.1) 
with quartic couplings fixed by supersymmetry at the scale Ms- The rg program 
uses now rg functions for the two Higgs doublet model to evolve the mass matrix 
for neutral scalars, plus the masses for charged Higgs bosons from Ms down to the 
electroweak scale. The mass of the lightest Higgs boson is smaller in this case. 

The fact that a clever choice of the renormalization scale permits to obtain 
reliable results in a leading log calculation can be used to derive simple analytical 
formulas for radiatively corrected Higgs masses (and couplings). This can be done 
by iteratively integrating rg equations to the required precision. 

4. Standard Model Stability Bound on the Higgs Mass 

In section 3 we saw that in the MSSM large radiative corrections to the lightest 
(CP even) Higgs mass arise in the case of a heavy supersymmetric scale. By using 
effective theory methods, the corrections were described as a renormalization group 
effect: the quartic Higgs coupling at the supersymmetric scale is small (as it is fixed 
by electroweak gauge couplings) but is driven to large values at the electroweak 
scale by the top quark contributions to (3\ ~ — 12/i^/(167r 2 ) which dominate the 
running of A if the top is heavy. 

This effect of top loop corrections is a purely Standard Model effect (the run- 
ning of A below the supersymmetric scale is described by the Standard Model rg 
functions) and has interesting consequences already in the pure Standard Model. 
The steepness in the running of A implies that if this quartic coupling is small at 
the electroweak scale (this means a light Higgs mass) it can be driven to negative 
values at a large scale if the top is sufficiently heavy. A negative value of A signals 
the appeareance of an instability in the effective potential at large values of the 
field: V((p) ~ A</? 4 — > — oo. If this pathology would appear for values of the field 
well beyond the Planck mass we shouldn't worry about it because we know already 
that the Standard Model has to be modified at such energy scales. In general, if the 
Standard Model is valid up to some large energy scale A (where some new physics 
will take over) we should be concerned about the possibility that the effective po- 
tential is destabilized below A. To avoid this pathology the Higgs mass should be 
heavy enough 

Mh > Mfj(M t ,A), 

where the critical value Mfj will be a monotonically increasing function of the top 
mass (the heavier the top is, the steeper the descent of A will be) and the scale 
A (the larger A is, the longer will be the running of A). The interesting point is 
that, for the current CDF/DO values of the top mass the limits Mfj are around 
100 GeV, a region with direct significance for future Higgs searches. It would then 
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be desirable to compute this critical masses, dubbed stability bounds, with good 
precision and this can be done following a procedure very similar to the one used 
in the previous section to compute the upper bound for the lightest MSSM Higgs 
boson mass. 

Before that, let us have a closer look at the instability of the potential. We 
need to compute reliably the value of the effective potential at large values of the 
field to see whether it is below the value at the electroweak minimum, in which case 
this minimum would get destabilized and eventually decay. As we don't know the 
exact potential we have to rely in perturbative approximations and so, a convenient 
choice of the renormalization scale should be made. To start just consider the choice 
[i ~ Mz that has the advantage that we know the values of the couplings without 
having to run them. The tree level potential is just 

with X(Mz) obviously positive, so that no instability would arise at high <p. If we 
add one-loop corrections we discover that the dominant piece comes from top loops 
[see fig. 7 (a)] and is 
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[Here rii counts the number of degrees of freedom of the i th particle with field- 
dependent mass M? = njip 2 (this form for the masses holds for the main contribu- 
tions to the potential. Higgs bosons have an additional field-independent piece, but 
their contribution is not important numerically) . Cj are some numerical constants.] 
We see that, for sufficiently large values of <p this piece dominates over the tree-level 
part and being negative drives the potential to negative values. Once again, the 
fact that AVi dominates over Vq does not imply that the perturbative expansion 
is not well behaved. In a similar way as what we discussed in the previous section 
the expansion parameter is of the form a log where again a ~ h 2 /16n 2 and now 
log ~ log(0//i). To be confident that the instability is really there, it has to appear 
for values of the field where a log < 1. The fact that this can actually happen is 
then a consequence of the hierarchy h t A of the relevant coupling in the loop 
correction to the relevant coupling in the tree level potential. On the other hand 
note that the log depends now on the scale \i. This means that fixing \i ~ Mz 
we cannot reliably study potential instabilities if they appear at field values much 
larger than (p ~ M z exp(47i). Assuming this is not the case, if a log < 1 but not too 
small one should care about higher order corrections. For example the dominant 
two-loop correction to the potential (at large (p) has the form 
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where the #'s are field independent constants (some combination of couplings) [the 
dominant diagrams are depicted in Fig. 7 (b,c)]. 

Exercise Taking advantage of the scale invariance of the effective potential compute 
the value of the constant # for the two-loop leading log contribution in (4.2) (keeping 
only h t and gs couplings). Use the knowledge of AVi and the one- loop rg functions 
leTrV 1 ) ~ Sh 2 and 16n 2 (3^ ~ h t [(9/2)h 2 - 8g 2 s ]. 




Fig. 7: Dominant loop diagrams contributing to the running of A. (a) One- loop 
and (b,c) two-loops . 

From the previous discussion it is clear that we can use our freedom in choosing 
ii to have a reliable perturbative expansion by the choice ii ~ ip. This choice will 
make small the radiative corrections so that we can just work with the tree level 
potential. Of course, this doesn't mean that the destabilization disappears (being 
caused by large loop corrections) . The size of radiative corrections is scale dependent 
and only the full potential doesn't depend on the scale. By choosing |i~^we just 
reorganize the expansion in such a way that the bulk of the effect is transferred from 
the higher loop corrections to the tree level part. The instability appears then as a 
result of A(/U ~ (p) < 0. We have then a clear picture of the procedure to follow: we 
consider the tree-level effective potential with parameters evaluated at some scale \i. 
When A(/u) runs to negative values, it means that the potential is being destabilized 
at <p of that order. Imposing that this destabilization scale is larger than the cut-off 
scale A we can compute the stability bound for that A: the lowest allowed value 
for A(/U ~ M z ) (and thus for M H ) is such that \{p, ~ A) = (which is in practice 
used as a boundary condition for the running A) . Implementing that program with 
one-loop rg functions amounts to compute the stability bound at leading-log order. 

For large A (say A ~ Mpi) the leading- log calculation can be a good enough 
approximation for the stability bound. A better precision can be achieved with a 
next-to-leading log calculation, proceeding along similar lines as what was done in 
section 3. Understandably, the NTLL calculation is mandatory in the case of lower 
A, say 1 — 10 TeV, when the effect of NTLL terms can become more important. 

The ingredients for a NTLL calculation should be clear by now. First, use a 
one-loop effective potential with parameters running with two-loop rg functions. 
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Next, keep control of the scale dependence of the potential. This has to be done 
in two regions of (p: in the region near the electroweak minimum, to ensure that 
the potential is minimized properly and the Higgs mass extracted correctly (as was 
done in the previous section) and then, in the region of the instability. In that 
region one can fix fj, = a<p and study the flatness of the bounds with the parameter 
a. Again, it turns out that, the tree-level calculation of the bound requires a careful 
choice of the scale while the one-loop calculation is much more stable. 

As a result of using a one-loop expression for the potential, one should replace 
the boundary condition A (A) = by a one-loop corrected version of it A (A) = 0. 
Here A is basically the one-loop coefficient of </? 4 (once \i = aip). From (4.1) 



log^-Ci 
or 



The fact that A should be used for a faithful description of the potential destabi- 
lization is exemplified by fig. 8 . For the indicated values of the parameters, the 
first plot gives the behaviour of the potential. Note the appeareance of a deep non- 
standard jmnimum at large values of the field. The second is a plot of the running 
of A and A with the scale, showing that the scale of destabilization of the potential 
indeed coincides with the scale at which A crosses zero. The point were A is zero is 
instead an order of magnitude lower. 




Fig. 8: Right: Effective potential. Left: Running A (dashed) and A (solid). 

Note that the potential does not run to — oo for large field values but rather 
turns again to positive values. This is caused by the decrease of ht at high scales. 
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The destabilizating force then weakens and eventually the effect of gauge couplings 
takes over making the couplings A and A positive again (see left plot in Fig. 8 ). 

The last ingredient for the NTLL computation of the bound is to extract phys- 
ical masses for Mh and M t as explained in Section 3 [Eqs. (3.10) and (3.11)]. 
The NTLL results are typically O(10) GeV lower than the LL ones and even more 
[O(20) GeV] for low cut-offs [e.g. A ~ 10 3 GeV}. The bounds, as a function of M t 
and A are plotted in fig. 9 . There is a non-negligible dependence on a s (larger for 
larger A) the bound being larger for smaller a s . 
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Fig. 9: Stability bounds for the Standard Model Higgs mass as a function of 
the top mass M t and the cut-off scale A from A = 10 3 GeV to 10 19 GeV in 
steps of two orders of magnitude. 

The area below a given curve corresponds to a region of parameters where the 
effective potential develops an instability at field values below the corresponding 
cut-off. This cut-off is usually referred to as the scale of new physics and should be 
interpreted with some care. The correct way of thinking about it is that, for a given 
pair (Mt, Mh) in the instability area, there is a curve labeled by some value of A 
that goes through that point. That value of A is the scale at which the instability 
develops in the pure Standard Model for the specified values of M t and Mh- As 
such, it can be computed very precisely, irrespective of whether A is large or small. 
In case that for the experimental value (M t , Mh) the instability scale A(M t , Mh) is 
below the Planck mass it would imply that the Standard Model would have to be 
extended in such a way that the instability is cured. This means that new physics 
should have an effect at scales below A. 
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5. Standard Model Metastability Bound on the Higgs Mass 

Strictly speaking, we may admit the possibility that the electroweak vacuum 
is not the deepest vacuum of the effective potential provided it is sufficiently long 
lived and the decay probability to the true vacuum is negligible. In the presence of 
a large potential barrier separating the two vacua it is reasonable to expect that the 
decay rate for tunneling (per unit time and unit volume) would be exponentially 
suppressed. However the Universe is old and large and it turns out that for some 
choices of the parameters below the stability curve ( e.g. for M t = 180 GeV and 
M/j = 100 GeV) the unstable electroweak vacuum should have decayed long ago 
and are thus unacceptable if no new physics is at work at (or below) the scale of 
the instability. Then, the requirement that an unstable vacuum is acceptable if its 
lifetime is longer than the age of the Universe weakens the stability bounds derived 
in the previous section increasing the allowed range of Higgs masses. 

However, there is a more stringent requirement an unstable electroweak vacuum 
should meet: it has to survive the high temperatures in the early Universe, when 
thermal excitations can trigger the decay to the true vacuum. The Higgs potential 
in the presence of the thermal plasma in the early Universe depends strongly on the 
temperature, T. For T much larger that the scale of instability, A, SU(2) x Z7 (1) 
symmetry is restored and the only minimum of the potential is at the origin. At 
T ~ A a new local minimum appears at a scale ~ T ~ A and at some temperature 
T c x the new minimum becomes degenerate with the one at the origin. For T < T c 
the minimum away from the origin is the deepest and the decay from the minimum 
at the origin to it becomes possible. For lower and lower temperatures the true 
minimum gets deeper and deeper evolving to the non-standard T = vacuum at 
T < A. The barrier between the two vacua is always present. Eventually, if the 
Universe remains trapped at the origin, at a temperature Tf? w of order 100 GeV 
the standard electroweak phase transition takes place, SU(2) x U(l) gets broken 
and the Universe sits in a new minimum that lies at the electroweak scale. The 
possibility of a no-return transition to the non-standard minimum is more likely than 
tunneling from quantum fluctuations at T = 0. For example, if M t = 180 GeV and 
Mh = 130 GeV, the unstable electroweak vacuum would be safe against quantum 
tunneling but cannot survive the high temperature early Universe. 

The phase transition to the non-standard minimum is strongly first order and 
proceeds via thermal nucleation of bubbles of true phase that grow till they fill the 
Universe. If a bubble of true phase is too small it collapses under surface tension 
while, if it is large enough the gain in potential energy compensates surface tension 
and the bubble grows. There is then a critical bubble, a saddle point of the free 
energy functional, that gives the critical energy that appears in the Boltzmann 
exponent for the decay rate. Consider a static and spherically symmetric bubble of 
true vacuum with Higgs profile <p(r). Here r is spatial distance from the center of the 
bubble, so that (f)(0) probes the instability region of the potential and (f)(r — > oo) 
goes to the false vacuum. At a given temperature, the extra energy for such a 
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configuration in a sea of false vacuum is given by 



E[d>(r)] = 4tt J™ r 2 dr \ {^j * + V(<f>, T) - V(vf a i se , T) 

The critical bubble minimizes this energy and so the critical bubble profile is the 
solution of the Lagrange equation 

d 2 <j> 2 # _ dV_ 
dr 2 r dr dfy 1 

with boundary conditions 

0(r -> oo) = v fa i se (T), d(j)/dr\ r=0 = 0. 

Here Vf a i se (T) is the field expectation value in the minimum close to the origin: 
Vfaise = for T > Tf w and w/ a z se = v E w(T) for T < Tf w . For a given tem- 
perature T, the effective potential is known and the critical bubble profile 0s(r) 
and energy Eb(T) = E[<f>B(r)] can be computed. The rate of nucleation of these 
bubbles per unit time and unit volume is then given by 

r/z/~wT 4 exp[-£ B (T)/T], (5-1) 

with uj a constant that can be taken of order unity for our purpose. The qualitative 
dependence of Eb with temperature is shown in fig. 10 . For Higgs masses larger 
than the LEP I experimental bound ~ 65 GeV the deepest minimum of the curve 
appears always in the region T 3> w so that the decay rate is maximal at 
temperatures much larger than Tf w . 

To get the probability of decay, the rate (5.1) should be multiplied by the 
volume of our current horizon scaled back to temperature T. The differential prob- 
ability is then 

dP M Pl 
-^ = K —e W (-E B /n 

with k ~ 3.25 x 10 86 and is plotted also in fig. 10 . Obviously it peaks near the 
minimum of E B /T. The integrated probability of decay is then 



I 



T - dP(T) 

~dT^ dT ■ (5 - 2) 



Note that this probability is not normalized to one. In fact its interpretation is that 
the fraction of space that remains in the old metastable phase is given by exp(-P). 
So, P ^> 1 (see fig. 11) indicates that the metastable electroweak vacuum would 
have decayed to the non-standard minimum in the hot early epoch of the Universe. 
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Fig. 10: Ratio of the critical bubble energy to the temperature (solid line) and 
differential probability for bubble nucleation (dashed) as a function of tempera- 
ture. T C EW and T\ are the two critical temperatures. 
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Fig. 11: Logarithm of the total probability of decay to the true vacuum as a 
function of Mh for Mt = 175 GeV . The point where the curve roughly crosses 
the dashed line gives the metastability bound. 
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Then, for fixed M t the probability P of thermal decay to the true vacuum can be 
computed as a function of M^. One example is shown in fig. 11 for M t = 175 GeV. 
Note the vast range of variation of P when is changed by a few GeV. Actually 
this is due to the exponential sensitivity of the decay rate to Eb and then to the 
shape of the effective potential. From fig. 11 we can derive the so called metastability 
bound for Mt = 175 GeV: below ~ 122 GeV P is very large and the vacuum decays 
quickly while for values above that, P is exponentially small and the electroweak 
vacuum, if metastable, is sufficiently long lived. The critical value P = 1 can be 
taken to compute the critical mass but as is clear from the figure the bound is 
insensitive to the exact value chosen for the critical probability [provided is very 
roughly 0(1)]. 

The metastability bounds would depend also on the cut-off scale A where new 
physics is expected to appear. By definition we can compute reliably the effective 
potential at field values lower than A and temperatures lower than Ta ~ A. The 
integrated probability that should be required to be less than 0(1) is then (5.2) 
with this temperature cut-off implemented 

P(T A )= [ A dP(T'). 
Jo 

For some fixed values of M t and M^, the condition P(T\) = 1 gives Ta. The profile 
of the critical bubble at that temperature can be calculated and in particular the 
value of the field at the center of the bubble </>s(0) obtained. To avoid the decay of 
the electroweak vacuum for this choice of parameters, new physics should modify 
the effective potential at values of the field of the order </>s(0), i.e. new physics 
should appear at A = <^b(0). 

The numerical results for the metastability bounds are shown in fig. 12 for 
two values of A: 10 19 GeV (three upper curves) and 10 TeV (three lower curves). 
For a given A the (M t , M^) plane gets divided in four different regions. Above the 
long dashed line the electroweak vacuum is absolutely stable, while below it is only 
metastable: there is a deeper non-standard minimum developing below A. Between 
the absolute stability lower bound and the solid line the metastable electroweak 
vacuum is long lived and acceptable. Below the solid line however, it would have 
decayed by thermal excitations in the early Universe. The short dashed line indi- 
cates the metastability bound for quantum tunneling at zero T. For A = 10 19 GeV 
the lower metastability bound on the Higgs mass is (for a s = 0.118) 

M h (GeV) > 2.306[M t (GeV) - 180] + 138. 

The implications of these bounds for the existence of new physics beyond the Stan- 
dard Model are discussed in the next and last section. 

t Again the precise meaning of this scale is provided by the calculation itself as the scale at 
which the effect of new physics should affect the finite temperature potential in a way suitable 
of modifying the decay rates. 



31 



l 1 1 1 1 1 1 1 1 1 r 




130 1 40 1 50 1 60 170 1 80 1 90 200 



M t (GeV) 

Fig. 12: Lower bounds on the Higgs mass for A = 10 19 GeV (upper set) 
and A = 10 TeV (lower set). Long dashed curves give the absolute stability 
bound, short dashed the metastability bound for quantum tunneling and solid 
for thermal tunneling in the hot early Universe. 



6. Implications 

The most important implication of the lower bounds on the Higgs mass derived 
in section 5 is that the measurement of Mh may provide an upper bound on the 
scale A of new physics. For this to be true, the top mass should be heavy enough 
so that some instability would appear in the Standard Model potential if the new 
physics were not present. Numerically the requirement is (all masses in GeV) 

M t >^ + 123GeV (6.1) 

(This is obtained from the metastability bound for A = 10 19 GeV). In particular, 
for Higgs masses above the experimental limit ~ 65 GeV the top quark should be 
heavier than ~ 152 GeV, which is the case. Similarly, (6.1) also tells that, unless 
M t < 160 GeV, the discovery of a Higgs boson by LEP II would imply that the 
Standard Model cannot be valid up to the Planck Scale. 

Moreover, in most of the parameter space, the scale below which new physics 
should enter is much smaller than Mp\. Consider for example fig. 13 where metasta- 
bility bounds for different top masses are plotted as a function of the scale A. Sup- 
pose now that LEP II finds a 90 GeV Higgs boson. Then, if M t > 170 GeV, we can 
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read from the figure that new physics should appear below A ~ 10 6 GeV. Other- 
wise, the Standard Model potential would be unstable and our vacuum would have 
decayed long ago. 
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Fig. 13: Metastability bounds on the SM Higgs boson mass as a function of the 
cut-off A for the indicated values of the top mass (lower curve for lower Mt). 
The horizontal line indicates roughly the reach of LEP II. 

Of course, the new physics should be such that this instability is cured. For 
example, a heavy fourth generation would make the effective potential even more 
unstable (a new large Yukawa coupling would add to the destabilization effect of 
the top quark) . Then, this kind of stability analysis can in fact be used to constrain 
such extensions of the Standard Model (see A. Novikov's contribution to these pro- 
ceedings). As a prime example of new physics that would go in the good direction 
we can consider the Minimal Supersymmetric Standard Model. In the MSSM, it 
turns out that the running of the quartic Higgs coupling is no longer dominated by 
the top Yukawa coupling. The effect of top quark loops is compensated by stops, 
the supersymmetric partners of the top. The relevant diagrams are shown in fig. 14 
. In fact, quartic Higgs couplings in the MSSM run like gauge couplings. This effect 
is schematically depicted in fig. 14 . Below the supersymmetric scale A the quartic 
Higgs coupling decreases steeply due to top loop effects. But once the supersym- 
metric threshold is crossed, all supersymmetric particles influence the running of A 
with the effect of Yukawas cancelling and leaving just the gauge renormalization. 
The figure then shows how the presence of the supersymmetric threshold prevents 
A to become negative and makes it rise again, thus stabilizing the potential. 
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In the framework of the MSSM as model for new physics, it is then tempting 
to try and put an upper bound on the supersymmetric scale by using metastability 
arguments. Unfortunately, as is shown already in fig. 13 , metastability bounds 
cannot compete with simple naturalness criteria that require Asusy < 0(1) TeV. 



Fig. 14: Schematic running of the quartic Higgs coupling below and above 
the supersymmetric threshold with the relevant diagramas indicated. The pair 
{Mf, Mfi) is chosen to lie in the instability region. 

Nevertheless, useful information can be extracted from the confrontation of the 
theoretical expectations for Higgs masses in the Standard Model and the MSSM. 
In fig. 15 we plot (steepest curve) the less restrictive metastability bound for the 
SM Higgs mass, choosing A = 10 19 GeV, i.e. assuming validity of the model up 
to the Planck scale. Superimposed in the same plot it is shown the absolute upper 
bound on the lightest Higgs mass in the MSSM for Asusy = 1 TeV (more or less 
the limit from naturalness). In both curves the uncertainty caused by as in the 
range indicated is reflected in the dashed curves. 

The (Mt, Mh) plane gets divided into four different zones. The experimental 
determination of the zone actually realized is of great interest. If the Higgs boson is 
found in the upper left region it would be compatible with the SM but it would be 
too heavy for the MSSM (remember that we are talking about a neutral Higgs boson 
with properties similar to the SM one). Alternatively, if it is the lower right area 
the one chosen by nature, then the Higgs mass would be too low for the pure SM 
and new physics should appear below the Planck mass. This new physics could well 
be in the form of the MSSM, and in that region the mass of the Higgs is compatible 
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with that hypothesis. Finally, the region in the lower left corner is compatible with 
the SM and MSSM while the upper right corner is incompatible with both. 




Fig. 15: Absolute lower met ast ability bound on the SM Higgs mass with A = 
10 19 GeV (steepest curves) and absolute upper bound on the lightest Higgs 
mass in the MSSM with Mg = 1 TeV and maximal stop mixing effect. Dashed 
lines show the variation when a s is changed in the indicated range. 

So, at least in some regions of parameter space, we may be able to distinguish 
indirectly between the lightest MSSM Higss boson in the decoupling SUSY limit 
and the SM Higgs, which is a very difficult task experimentally, fn any case it is 
remarkable that the detection of a Higgs boson, the last missing ingredient of the 
Standard Model, can already provide information against the pure SM and point 
towards the need of new physics at some scale well below the Planck scale. 
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